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Abstract 

A new approach developed in Ref.l for investigating the spatial boson corre- 
lations at low temperatures is considered in the particular case of the Coulomb 
interaction potential. This approach is based on integro-differential equations 
for the radial distribution function and can be used beyond the weak coupling 
regime. 

PACS: 05.30.Jp, 67.90. +z 

Description of the charged Bose liquid, or, in other words, of the strongly correlated 
many-boson system with the Coulomb interparticle potential, is a long-standing problem 
of the condensed matter physics. Recently, interest in this problem has been renewed in the 
connection with the noteworthy results of applying the Bose and Bose-Fermi liquid mod- 
els [§] to the explanation of various phenomena inherent in the superconductors with high 
critical temperatures. These models are based on the Bogoliubov approach || implying the 
weak coupling regime. For the system of charged bosons the weak coupling approximation 
can be used at small values of the Brueckner parameter r# ~ 0.1 . However, questions re- 
quiring detailed investigations often concern the materials with relatively low densities of 
charge carriers which correspond to r^— values even larger than those of the electron gas at 



the metallic densities. For example, there are the results |4j] testifying to the possibility of 
forming the so-called large bipolaron (e.g. the two-electron cluster with a large radius of 
electron-electron separation, see papers in Ref. |J) inside the polaron environment at the 
electron densities ~ 10 19 cm -3 close to the carrier densities of the heavily doped degener- 
ate polar semiconductors. To investigate the possible phase transition of the charge carrier 
system from the polaron state to the bipolaron phase at these densities, one should operate 
with the charged Bose liquid of large bipolarons for r$ ^> 1 . Indeed, the " large-bipolaron " 
Bohr radius can be expressed f§] as 
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where Eq is the static dielectric constant of the polar material considered, mu v denotes the 
bipolaron mass, and e is the electron charge. To connect au p with the canonical electron 
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Bohr radius a e = = 5.29 • 10 cm , we use the reasonable estimate H mu v ~ 2 m po \ 
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(m po i is the polaron mass) and the relation m po i ~ — — - m* (m* is the band mass) 

correct at the values of the Frohlich electron-phonon coupling constant a ~ 7 typical for 

Ti 2 

the large bipolaron H. This leads to au p ~ 0.0026 e r . Further, with the tables of the 

physical properties of polar materials |J we find that Eq ~ 10 4- 100 and m* e ~ 0.1-r 0.5 m e 

represent the most frequent case. Choosing, say, Eo = 100 and m* = 0.1 m e , we obtain 

cibip ~ 2.6a e = 1.375 ■ 10 -8 cm . With this choice the "large-bipolaron" Brueckner parameter 
1 ( 3 V /3 u i i. i i 19 3 



r$ = has values ~ 20 at the bipolaron densities n ~ 10 cm , which proves 
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the statement mentioned above. Note that another choice gives even larger values of r§ . 
Thus, to make the final conclusions on the possibility of the phase transition from the Fermi 
state to the Bose one in the charge carrier systems of various materials and to clear up the 
status of the Bose liquid models [Q], one should deal with the Bose many-particle systems 
beyond the weak coupling regime. 

The aim of this letter is to consider potentialities of a new approach to investigating 
the spatial correlations of cold bosons in the particular case of the Coulomb interaction 
potential. The approach is based on integro-differential equations for the function of the 



radial boson distribution that have been derived in paper Q in a way analogous to that of 
Ref. J7|. For the charged bosons the simplest of these equations is written as 

—^prAg^ir) = ^ + nf yS^Uy) ~ l)dff. (1) 
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where g(r) denotes the radial distribution function, m and Q are the boson mass and charge, 
n is the particle density. Besides, A stands for the Laplacian. To be convinced of the capacity 
of the equation ([I]) for giving reasonable estimates of the radial distribution function of the 
charged bosons, let us consider the weak coupling regime characterized by 

g(r) = 1 -e(r), e(r) < 1. (2) 

Using (0) and only keeping the terms linear in e(r) , we arrive at the linear version of the 
equation (JTJ) : 
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which can easily be solved via the Fourier transformation. The Fourier transform of e(r) 
obeys the relation 
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Introducing the new variable x = q/A with A 4 = IQirnmQ 2 / h 2 , from (^) we find 

() = ^(x) , x 4 
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It is now interesting to compare the estimate (|5]) with the result M of the random phase 
approximation (RPA) 

-/ x Irpa{x) x 2 
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Short-range correlations are determined by the behaviour of e(q) at x = q/A ^> 1 for which 
we have 

fest{x) ^ , /kPa(z) ~ ^-j (X > 1). 



There is also good agreement between f es t(x) and f RPA (x) in the region x ~ 1 ruling the 
spatial correlations at the intermediate distances close to the screening radius 1/A: 

f est (l.5) = 0.090 , f RPA (1.5) = 0.086 ; 

/ est (l) = 0.333, f RPA (0.5) = 0.293; 
/ est (0.5) = 0.889 , f RPA (0.5) = 0.758 ; 
When i < 1 (long-range correlations), we again obtain close numerical results 

f est {0.2) = 0.997 , f RPA (0.2) = 0.960 , 

however, f es t(x) and f RPA (x) tend to unity in completely different ways: 

f est (x) ~ 1 - 2 x 4 , f RPA (x) ~l-x 2 {x < 1) . 

The comparative consideration of f es t{x) and f RPA (x) allows one to make the following two 
conclusions. The first, negative, is that the model presented is not able to yield the correct 
behaviour of the structure factor S(k) = 1 — ne(k) at small k and, therefore, the correct 
spectrum of the low-lying excitations given by k 2 / S(k) (see Ref. [[]]) . The second conclusion, 
positive, consists in that the equation (JIJ) can provide a good estimate of the pair correlation 
function h(r) — g(r) — 1 for distances r C 1/A and r ~ 1/A exactly where the spatial boson 
correlations are significant. 

These conclusions have been found in the weak coupling regime. And what can we 
expect, as to the results of ([]]), beyond the weak coupling, for r s > 1. It is known |10 



that RPA should be corrected for rs > 1 to avoid the trouble of negative values of g(r) at 
small particle separations. This occurs because the short-range correlations are taken into 
account not in the best way within RPA. At small enough the short-range correlations 
do not yield any essential contribution to the thermodynamic quantities. Therefore, this 
shortcoming of RPA does not practically affect the results for rs < 1 ■ On the contrary, at 
large rs the short-range correlations become significant and influence the thermodynamic 
quantities considerably. In this case RPA leads not only to the negative values of g(r) at 



small r but overestimates the absolute value of the ground state energy. As to the model 
based on the equation (JI|), it has no problem like this. Indeed, at small r where 
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the expression (JT[) can be written as 

--A^ 2 (r) + ^^ 2 (r) = 0. (7) 
m r 

It is nothing else but the Schrodinger equation like those which are often used to investigate 
the spatial boson correlations at small particle separations as well as in dilute systems [|TJ. 
Thus, we may expect that the relation (P will also give reasonable results on the charged 
boson correlations for rg > 1. 

If we have a reliable estimate of g(r), we can obtain adequate estimates of the thermo- 
dynamic quantities of the system considered. In particular, it is possible to calculate E , the 
ground state energy of the charged bosons in the neutralizing background, with the following 



well-known relation [11 





Here N denotes the number of bosons, is the energy of noninteracting particles (in our 
case Eid = ). The quantity 7 is the coupling constant and g(r, 7) denotes the solution 
of the equation obtained with the replacement Q 2 — > jQ 2 in ([!]). On the basis of the 
relation (||) we find 

which can be rewritten as 
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Substituting (0) into ® we obtain 
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The result of numerical integration 
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J dl J dy ysm(ly) 1- 4 Q J = 0.9339 
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makes it possible to derive from (|TTJ) the relation 

E 



= -0.2378 Q 2 A = -0.885 r 5 3/4 Ry (12) 



where rs and Ry obey the standard definitions 

mQ 2 ( 3 \ 1/3 „ mQ 4 

The RPA result for E/N can easily be found with the relation (|6]). This yields the expression 
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Numerical integration 



J d 1 j dl J dy ^ysm(ly) (l - j . (13) 
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J dl J dy ysm(ly) h - _|L= j = 0.8468 
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allows one to obtain 



77 

— = -0.2156 Q 2 A = -0.803 r~ 3/4 Ry . (14) 



Thus, agreement between the results of (0) and fll4]) is within the accuracy of 10% . Besides, 
it is also interesting to make a remark on the next-to-leading contribution to the mean 
energy of the system. Within the approach considered this contribution can be calculated 
substituting 

g (r) = 1 - £i(r) - e a (r) (e 2 (r) < e x (r) < 1 ) 
into the equation (JJ). Omitting details of the calculations, let us only write the final result 
E 

— = (-0.885 r s 3/i + const) Ry , canst w 0.08 (15) 
which can be compared with the relation 



_ = (-0.803 r s 3/4 + 0.212) Ry (16) 

found within the Bogoliubov scheme 0. 

Summarizing, let us recall the basic points of the letter once more. The new approach 
of investigating the spatial correlations of cold bosons has been considered in the particular 
situation of a many-boson system with the Coulomb interaction potential. In this case the 
radial distribution function g(r) can be estimated with the integro-differential equation (ffj). 
As it has been shown, this equation is able to yield reasonable results on the short- and 
intermediate-range correlations. However, the obtained decay of the function h(r) = g(r) — l 
at distances essentially larger than the screening radius is not correct. Nevertheless, Eq.(|l]) 
has provided a good estimate of the ground state energy at small r$ . This result gives an 
optimistic view on the possibility of using Eq.(|l]) to explore the thermodynamic behavior of 
the charged Bose liquid ( r$ > 1 ) where the influence of the short-range correlations on the 
system properties is essentially more significant than that for rs < 1 ■ 
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